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It is well known that a naked singularity occurs in the gravitational collapse of an inhomogeneous 
dust ball from an initial density profile which is physically reasonable. In this paper we show that 
explosive radiation is emitted during the formation process of the naked singularity while we fix the 
background spacetime. The energy flux is proportional to (tcH — i) -3 ^ 2 for a minimally coupled 
massless scalar field, while it is proportional to (tcH f° r a conformally coupled massless scalar 

field, where ten — t is the "remaineing time" until the distant observer could observe the singularity 
if the naked singularity was formed. As a consequence, the radiated energy grows unboundedly 
for both scalar fields. The amount of the power and energy depends on the parameters which 
characterize the initial density profile but do not depend on the gravitational mass of the cloud. In 
particular, there is a characteristic frequency v s of the singularity above which the divergent energy 
is radiated. The energy flux is dominated by particles of which the wavelength is about ten — t at 
each moment. The observed total spectrum is nonthermal, i.e., vdN/du ~ (v/v a )~ l for v > v a . If 
the naked singularity formation could continue until a considerable fraction of the total energy of the 
dust cloud is radiated, the radiated energy would reach about 1O 54 (M/M0)erg. The calculations are 
based on the geometrical optics approximation which turns out to be consistent for a rough order 
estimate. The analysis does not depend on whether or not the naked singularity occurs in its exact 
meaning. This phenomenon may provide a new candidate for a source of ultrahigh energy cosmic 
rays or a central engine of 7-ray bursts. 

PACS numbers: 04.20.Dw, 04.70.Dy, 98.70.Sa 



I. INTRODUCTION 



It was shown that the collapse of an inhomogeneous dust ball results in shell-focusing naked singularity forma- 
tion |l| . Ori and Piran and Harada showed that naked singularity formation is possible from the spherical 
collapse of a perfect fluid with a very soft equation of state. Moreover, a kind of runaway collapse in Newtonian 
gravity looks like naked singularity formation in general relativity. These strongly suggest that the collapse in the 
Universe will lead to a drastic growth of spacetime curvature which is not covered by the event horizon. 

In this respect, a number of researchers have examined the emission during naked singularity formation. In classical 
theory, Shapiro and Teukolsky |5|,[l0| reported the drastic growth of a curvature polynomial but no apparent horizon 
in their numerical simulation of the axisymmetric collapse of a prolate spheroid of collisionlcss particles. They also 
reported that little gravitational radiation was emitted in their simulation. On the other hand, Nakamura, Shibata 
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and Nakao mU suggested that the naked singularity formation in the collapse of a prolate spheroid may provide a 



strong source of gravitational waves. Recently, Iguchi, Nakao, and Harada |12|] and Iguchi, Harada, and Nakao 113 14 



examined the behavior of nonspherical linear perturbations of the collapse of an inhomogeneous dust ball. One of 
their main results is that perturbations of a curvature polynomial grow unboundedly while the energy flux of the 
gravitational waves remains finite as the Cauchy horizon is approached, which implies that the instability is rather 
milder. 

On the other hand, Hawking showed that thermal radiation is emitted in the gravitational collapse to a black 
hole. The emission can be interpreted as particle creation due to the time variation of a strong spacetime curvature. 
In the collapse to a black hole, the magnitude of the spacetime curvature of the observable region is, at most, the 
order of the inverse square of its gravitational mass because the singularity is covered by the event horizon. On 
the other hand, in the collapse to naked singularity, the magnitude of the spacetime curvature outside the event 
horizon grows unboundedly and the strongly curved region can be seen by an observer. This fact suggests that 
naked singularity formation may be a strong source of emission due to the quantum effect. In this context, Ford 
and Parker Jl6| calculated the radiation during the formation process of shell-crossing naked singularity. They found 
that both the power and the energy remain finite, which is because of the weakness of the shell-crossing singularity. 
They also showed that the collapse of supercritically charged matter to the naked singularity emitts unbounded power 
and energy, although this spacetime is not a solution of the Einstein equation. Hiscock, Williams and Eardley JTt]] 
showed the diverging energy flux in the self-similar collapse of a null dust and also calculated the expectation value 
of stress-energy tensor in its two-dimensional version. Barve et al. []l8f showed the diverging power of radiation in 
the self-similar collapse of a dust ball. They also calculated the expectation value of stress-energy tensor in its two- 
dimensional version [19| . Vaz and Witten |2(]] also calculated the Bogoliubov coefficients which relate to the spectrum 
of emission for this model. The above studies were all based on the geometrical optics approximation. 

The last two examples in which the diverging power is emitted are self-similar collapse. However, the self-similar 
collapse is a particular solution among gravitational collapse solutions in general relativity. Moreover, we do not 
know whether or not the central part of the realistic collapse of the nonzero pressure perfect fluid tends to be self 
similar in a strong-gravity regime such as a shell-focusing singularity, although such tendency has been observed in 
the Newtonian regime. For the self-similar collapse of a dust ball, it has been shown that the redshift at the center 
becomes infinite and that the curvature strength of the naked singularity is very strong. Moreover, it is known that 
such a solution does not allow the initial density profile which is a C°° function with respect to the local Cartesian 
coordinates. On the other hand, Christodoulou |2j imposed the C°° condition on the initial density distribution in 
his famous paper on the violation of cosmic censorship in the collapse of a dust ball. For the C°° case, the redshift 
is finite and the naked singularity is not very strong. See ]2l|-|23[| for details. Though the Einstein equation does not 
require such strong differentiability to initial data, we usually set such initial data in most astrophysical numerical 
simulations. Although this model is shown to be unstable to nonspherical perturbations, the model is still good for 
the sufficiently spherically symmetric collapse of dust until the nonspherical perturbations grow sufficiently. It will 
be important to examine whether or not the divergence of the radiated energy flux in the self-similar collapse only 
comes from the uniqueness of self-similarity. In this context, the authors have already reported explosive radiation 
during the formation process of the naked singularity which occurs in the collapse of an inhomogeneous dust ball from 
the initial density profile which is a C°° function [|24j| . In this paper, we present the details of the calculations in the 
above paper and examine power, energy, and spectrum of that radiation process. 

This paper is organized as follows. In Sec. [Il|, we will recall the least of quantum particle creation by a collapsing 
ball based on the geometrical optics approximation. In Sec. [II, we present the background solution of spacetime 
which describes the collapse of an inhomogeneous dust ball. In Sec. IV, we present numerical results of power and 
energy due to the quantum particle creation and give some physical discussions. In Sec. [v|, we determine the spectrum 
of the radiation which is important not only astrophysically but also in estimating the validity of the geometrical 
optics approximation. In Sec. VI, we analytically examine n ull g eodesies and redshifts in this spacetime and thereby 
understand some featu res of the numerical results. In Sec. VII, we estimate the validity of the geometrical optics 
approximation. In Sec. VIII , some implications of the emission from the forming naked singularity are discussed. In 
Sec. IX , we conclude the paper. We use the units in which G = c = % = 1 unless explicitly stated. We follow the sign 
conventions of Misner, Thorne, and Wheeler [E5|. 



II. PARTICLE CREATION BY A COLLAPSING BALL 



We consider both minimally and conformally coupled massless scalar fields in a four-dimensional spacetime which is 
spherically symmetric and asymptotically flat. Let T, R, 6, <fi denote the usual quasi-Minkowskian time and spherical 
coordinates, which are asymptotically related to null coordinates u and v by u ps T — R and v ~ T + R. An incoming 
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null ray v = const, originating on propagates through the geometry becoming an outgoing null ray u = const, 
and arriving on S + at a value u = F(y). Conversely, we can trace a null ray from u on 3 + to v = G(u) on 
where G is the inverse of F. Here, we assume that the geometrical optics approximation is valid. The validity of 
this approximation will be discussed later. The geometrical optics approximation implies that the trajectories of the 
null rays give a surface of the constant phase. Then, in the asymptotic region, the mode function which contains an 
ingoing mode of the standard form on 3~ is the following: 
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(2.1) 



where we have imposed the reflection symmetry condition at the center. On the other hand, in the asymptotic region, 
the mode function which contains an outgoing mode of the standard form on ^ + is the following: 
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Note that in the above we have normalized the mode functions as 

(u w i m ,u u q> m ') = 8(uj - uj')8 w 8 mm >, 
where the inner product is defined by integration on the spacelike hypersurface S as 



Using the above mode functions we can express the scalar field (f> as 
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According to the usual procedure of canonical quantization, we obtain the following commutation relations 

5{lo - uj')5w5 mm >, 
5(uj - uj')5w5 mm >, 
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(2.3) 
(2.4) 

(2.5) 
(2.6) 

(2.7) 
(2.8) 



where it is noted that the Lagrangian in the Minkowski spacetime is common for both minimally and conformally 
coupled scalar fields. Here, a™ lm and a™^ are interpreted as annihilation operators corresponding to in and out 
modes, respectively. Then we set the initial quantum state to in vacuum, i.e., 



<CJo >= o. 



(2.9) 



The power for fixed I and m is given by estimating the expectation value of the stress-energy tensor through the 
point-splitting regularization in a flat spacetime as U& 
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for a minimally coupled scalar field, and 
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(2.10) 
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for a conformally coupled scalar field. It implies that the power depends on the method of coupling of the scalar field 
with gravity. However, if 



G 7 



~G> 



(2.12) 
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holds, the radiated energy from u = atow = &ofa minimally coupled field 



and that of a conformally coupled field 



rb 

E lm = I P lm du (2.13) 



E lm = I P lm du (2.14) 



coincide exactly. The total power is given by summation of all {l,m). The simple summation diverges. This is 
because we have neglected the back scattering effect by the curvature potential, which will reduce the ra diated flux 



consid erably for larger I. Therefore we should recognize that the above expressions of the power, Eqs. ( 2.1C ) and 
( |2.11 ), are a good approximation only for smaller I. Hereafter we omit the suffix I and m. 



The spectrum of radiation is derived from the Bogoliubov coefficients which relate the in and out modes given 

by© 
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From the above equations, we find the following useful relation: 

Pw'w = —ia u >(-w)- ( 2 - 17 ) 
The Bogoliubov coefficients satisfy the following relations: 
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dCu{a G]ll] al UJ , - = <K W - ( 2 - 18 ) 

d£j{ac }u} Ps>u> - PuwOiQu') = 0. (2-19) 
The expectation value N(u>) of the particle number of frequency lo on 3 + is obtained by 

/>oo 

N(u) = / duj'\^\ 2 . (2.20) 



It is noted that these results are free of ambiguity coming from local curvature because the regularization is done 
in a flat spacetime. 



III. LEMAiTRE-TOLMAN-BONDI SOLUTION 

The spherically symmetric collapse of a dust fluid is exactly solved (2^j2^|. The solution is called the Lemaitre- 
Tolman-Bondi (LTB) solution. The metric is given by 

ds 2 = -dt 2 + fo^k^ dr 2 + R 2 (t,r)(d0 2 + sin 2 9dtf) (3.1) 

in the synchronous comoving coordinates. The matter field is given by 

T^ v = eu»u\ (3.2) 
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F(r) and f(r) are arbitrary functions. F(r) is the mass function, which is twice the Misner-Sharp mass. f(r) is the 
energy function. The areal radius R satisfies the following equation: 
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(R,) 2 = ^ + f(r). 



Equation (3.4) is integrated to the following form 

t-t.(r) = 

where g(y) is a positive function given by 
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for < y < 1, 
for y = 0, 
for y < 0, 



(3.6) 



t s (r) is a constant of integration, and the sign is chosen so that the collapsing phase is relevant. Note that the metric 
has physical meaning only on < r and t < t s (r). We rescale the radial coordinate r as follows 



R(0,r) = r. 



(3.7) 



From Eq. ( |3.3| ), if we require that the initial density profile at t = is a C°° function with respect to the local 
Cartesian coordinates, F(r) should be expanded around r = as 



F(r) = F 3 r 3 + F 5 r 5 + F 7 r 7 + ■■■. 
where we assume that F3 is positive. Then, t s (r) is given by 
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We can easily find that the time of the occurrence of shell- focusing singularity is given by t s (r), while the time of the 
occurrence of apparent horizon is given by 



tah(r)=t s (r)-Fg(-f). 

Hereafter we assume marginally bound collapse (/ = 0) for simplicity. 

At an arbitrary radius r, the LTB spacetime can be matched with the Schwarzschild spacetime 
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The matching conditions with the Schwarzschild spacetime are obtained by requiring that the surface moves on the 
timelike geodesic and that the three metric on the timelike hypersurface should be continuous. The result is 
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(3.12) 
(3.13) 
(3.14) 



at the surface r = r s r. The retarded time u and the advanced time v in the Schwarzschild spacetime are defined by 



u = T-R- 2M In 
v = T + R + 2M\n 



R - 2M 

2M ' 
R-2M 

2M 



(3.15) 
(3.16) 



For the case of marginally bound collapse, R is given by 
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(3.17) 



The occurrence of singularity is at the time 

ts(r) = -^J. (3.18) 

Each mass shell is trapped at the time 




(3.19) 



We denote the time of the occurrence of the singularity at the center as to = t s (0). 

Detailed analysis on the occurrence of naked shell-focusing singularity has been done in (l]-^|,|2^ -^J. Here we 
pres ent one of the results. If we assume marginally bound collapse and assume the initial density profile expressed by 
Eq. ([0]), then the central shell- focusing singularity is naked if and only if F5 is negative. The singularity satisfies not 
the strong curvature condition (SCC) defined by Tipler |}2| but only the limiting focusing condition (LFC) defined 
by Krolak []33"f for radial null geodesies which terminate at the singularity. On the other hand, the singularity satisfies 
both LFC and SCC for radial timelike geodesies pjj . 



IV. POWER AND ENERGY 



A. Globally naked singularity 

We can find the function G or F by solving the trajectories of outgoing a nd in going null rays in the dust cloud 



and determining the retarded time u and the advanced time v through Eqs. ( 3.15 ) and ( 3.16 ) at the time when the 
outgoing and ingoing null rays reach the surface boundary, respectively. The trajectories of null rays in the dust cloud 
are given by the following ordinary differential equation: 

^ = ±R,r, (4.1) 
ar 

where the upper and lower signs denote outgoing an d in going null rays, respectively. 

We have solved the ordinary differential equation (fO]) numerically. We have used the Runge-Kutta method of the 
fourth order. We have executed the quadruple precision calculations. 

We choose the mass function as 

F{r) = F 3 r 3 +F 5 r 5 . (4.2) 

We find that the central singularity is globally naked for very small r s f if we fix the values of F3 and F§. Although 
we have calculated several models, we only display the numerical results for the model with F3 = 1, F5 = —2 and 
r s f = 0.02 in an arbitrary unit because the features are the same if there exists globally naked singularity in the 
model. The total gravitational mass M is given by M = 3.9968 x 10~ 6 for this model. See Fig. |l]for trajectories of 
null geodesies. We also indicate the location of the singularity, the apparent horizon, and the Cauchy horizon in this 
figure. 

See Fig. ^| In Fig. ^|(a), the relation between u and v, i.e., the function G(u) is shown, where uq (vq) is defined as the 
retarded (advanced) time of the earliest light ray which originates from (terminates at) the singularity. In Fig. ||(b), 
the first derivative G'(u) is shown. This implies that G'(u) does not diverge but converges to some positive constant 
A with < A < 1. In Fig. ||(c), it is found that the second derivative G"(u) does diverge as u — > uq. This figure 
shows that the behaviors of the growth of G"(u) are different from each other during early times (10 -4 < Uq — u) and 
during late times (0 < uq — u ^ 10~ 4 ). During early times, the dependence on u is written as 

G" cx -(-u)- 2 , (4.3) 

while, during late times, the dependence is written as 

G" cx -{uq - u)~ 1/2 . (4.4) 
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Here we determine the magnitude of the power by physical discussions. We assume that the particle creation during 
early times is due to the collapse of the dust cloud as a whole while that during late times is due to the growth of 
the central curvature. Then, we consider a constant of proportion which will appear in Eqs. (4.3) and (4.4). First 
we should note that the coefficient must be written by initial data because it must not depend on time. Next we 
should note that we can identify any t — const < to hypcrsurface with the initial hypersurface. Since the constant 
of proportion directly enters the expression of the power which is phy sically me aningful, the coefficient must be 
independent of the choice of an initial slice. On the other hand, Eqs. (4.3) and (4.4) demand a dimcnsionful constant 
of proportion. For early times, the only possible quantity is the gravitational mass M of the dust cloud. On the other 
hand, for late times, as derived in Appendix [A| the only possible quantity is (ioAo); where to is given by 



to - 3F3 



1/2 



and Iq denotes the scale of inhomogeneity defined as 



-F 5 



-1/2 



Thus we can determine the coefficient except for a numerical factor as 

G" w -f e M(-u)- 2 , 
for early times, while the dependence is written as 

-1/2 



G" « -fiA 



(uo-u)- 1 / 2 , 



(4.5) 



(4.6) 



(4.7) 



(4.8) 



for late times, where f e and fi are dimensionless positive constants of order unity. It implies that there is a charac- 
teristic frequency of singularity which is defined as 



(4.9) 



The early time behavior derived above is good for — u 3> M, while the late time behavior is good for < uq~ u <C uj s . 
The turning point from the early time behavior to the late time behavior is roughly estimated as 



) 2/3 ^r 1 



(Muj s )- 1/3 M. 



Uq — u « {Mlj s 
Here we define the turning point frequency u)t P as 

uj tp = (Mu s )- 2 '*u a = {MtOsY^Nr 



(4.10) 



(4.11) 



The above estimates show a good agreement with the numerical results. 

Since G" and therefore G'" diverge as u — > uq, the power of radiation diverges for both minimally and confo rmall y 
coupled scalar fields. For this model, we have numerically calculated the power of radiation by Eqs. ( 2.1 0| ) and (2.11). 
The results are displayed in Fig. ||. 

Ba sed o n the abov e physical discussions, we can obtain the formula for the power by the particle creation using 
Eqs. ( |2.10D and (|2.1l| ). As seen in Fig. |(b), we find 



G' w 1 

during early times, and 

G' m A, 

during late times. Then, the power during early times is obtained as 

P - ^/eMC-u)-, 



(4.12) 
(4.13) 

(4.14) 
(4.15) 
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The power during late times is obtained as 



P^^f^iuo-uy^, (4.16) 
P « - a) -1 - (4-17) 

Therefore, the power diverges to positive infinity for both minimally and conformally coupled scalar fields as u — ► ito- 
The radiated energy is estimated by integrating the power with time. During late times, the radiated energy is 
estimated as 

E " 2S / " J . 1/2 <"°-">~" 2 - <4 - 1S) 
E^L^^l (4,9) 

Therefore, the total radiated energy diverges to positive infinity for both minimally and conformally coupled scalar 
fields as u — > Uq. However, in realistic situations, we may assume that the naked singularity formation is prevented 
by some mechanism and that the quantum particle creation is ceased at the time uq — uk At. In other words, G" /G' 
tends to vanish for uq — u < At. In such situations, the total divergence term in the expression of the power of a 
minimally coupled scalar field gives no contribution to the total radiated energy. Therefore, the total energy for a 
minimally coupled scalar field and that for a conformally coupled one coincide exactly, i.e., 

The relation of At and the maximum central energy density ^ ax of the dust cloud at that time is given by 

e rnax _ ( A £)- 2 . (4.21) 

From the above numerical results and physical discussions, we obtain the following formula for the function G(u) 
or F(v): 

G(u) w u + f e M\a(-u) (4.22) 
F(v) re v - f e M ln(-u) (4.23) 



or 



for early times, and 



G(u) « A(u - uo) - -Af^iuo - uf/ 2 + v (4.24) 



F(v) « A-\v - v ) + ^fiojV^A-Hvo ~ v)Y' 2 + u (4.25) 



for late times. 



B. Black hole 



For comparison and for a test of the numerical code, we have also calculated the function G for the Op penh cimcr- 
Snyder collapse. The model is given by setting F$ = 1, F$ = 0, and r s t — 0.2 (M = 4 x 10 -3 ) in Eq. (4.2) in an 
arbitrary unit. 

The final fate of the Oppenheimer-Snyder collapse is a spacelike singularity covered by the event horizon. Then it is 
expected that radiation during the collapse tends to Hawking radiation as u — > oo. In reality, this is true. See Fig. ^ 
in which trajectories of null geodesies are displayed. We can see that the event horizon covers the singularity. See 
Fig. H for the behavior of the functions G(u), G'(u), and G"(u). During early times, the behavior of G(u) is basically 
the same as that of the globally naked singularity, i.e., G(u) is written as 
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G(u) ku + f e M ln(-w). (4.26) 

This is very reasonable because the behavior during early times will be independent of details of the density profile 
within the dust cloud. During late times, G is written as follows 

G(u) w -const x exp f j^J + v h , (4.27) 

where an ingoing null ray with Vh is reflected to the outgoing null ray which is on the event horizon. This relation 
will be derived later. 

The power of radiation has been calculated numerically. During early times, the pow er is the same as that for 
globally naked singularity. The results are shown in Fig. g. On the other hand, from Eq. (4.27), we can estimate the 
power and the radiated energy during late times as 

The numerical results show a good agreement with these analytic results. The numerical result for the power of 



radiation in this model has reproduced that of the Hawking radiation (4.2S) by a black hole within ~ 0.01 % accuracy. 

We should also note that the late time radiation during the LTB collapse with the locally naked singularity is the 
same as that of the Oppenheimer-Snyder collapse. This is because a distant observer cannot see the locally naked 
singularity. 



V. SPECTRUM 



A. Total spectrum 



Now that we have obtained the function F(v), we can calculate the spectrum of radiation by Eqs. ( 2.15 ), fl2.16 ), and 
(2.20). Some analytic expressions for the Bogoliubov coefficients are obtained in Appendix |b| In order to determine 
the spectrum nume ricall y, we i ntrod uce the Gaussian window function in the Fourier transformation as usual. That 
is, in place of Eqs. fl2.15| ) and ( |2.16| ), we have calculated the Bogoliubov coefficients 
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F(v)- 



exp 



F(v) - u c 



(5.1) 
(5.2) 



where u c and £ are the center and the width of the Gaussian window function, respectively. Since the wave packet 

2" 
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(5.3) 



has the frequency band width Auj ~ 2/S, we can obtain the number of particles per unit frequency band as 



dN 



2 



f °° duj'lP^ 



(5.4) 



where we have again assumed that the geometrical optics approximation is valid. 

Although we could use the calculated data for the function F(v), we have used the analytic formula which have been 
derived based on the numerical results for convenience of retaining accuracy. Since the diverging power is associated 
with the late time behavior, we concentrate on the late time radiation. In order to obtain the total spectrum, wc 
extrapolate the function F linearly as 
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(A x (v - v ) + u a (v < v), 

F{V)= ) A-i(v-v ) + ^f l ^{A-\v -v)}^ + u a ( Vl <v<v ), ^ 



where V\ is given by 



3" 



V - Vi 



(i-^)zr 1 " 



Awj 1 . (5-6) 



We should note that Eq. ( |5.5| ) can be used for u — u < lo^ p . Later we will show that Lo tp ^ w s for r s / <C Iq. In fact, 
it turns out that we only have to pay attention to the spectrum above uj tp . We have chosen u c — u and £ to be 
rather smaller than u — Ux with u\ = F(y\), for the contribution dominantly comes from u f=a Uq. We should note 
that there is no radiation for u > uq in this extrapolation. We could adopt another extrapolation, for example, an 
extrapolation such that F(v) — u$ is antisymmetric with respect to v = vq, i.e., 

, A~ x {v - v ) - -fiul'^A- 1 ^ - w )] 3/2 + u {v <v< 2v - vi), 
F(v) = { | (5.7) 

A~ x (v - v ) + -fiujl^A-^vo - v)} 3 / 2 +u a ( Vl <v< w ), 

This a ntisy mmetric extrapolation turns out to only double the amplitude of the spectrum through the linear extrapo- 



lation (5.5). Therefore, we only present the spectrum through the linear extrapolation. Because the second derivative 
of F is diverging as v — > vo, we can only require that the first derivative of F should be continuous at v = vq. If 
we allow discontinuity of the first derivative, radiation due to this discontinuity dominates the radiated energy flux, 
which is out of concern in this paper. 

The obtained spectrum is shown in Fig. |^. The parameters are fixed as A — 0.8 and fi = 1. We should note that 
the contribution to the total radiated energy mainly comes from lo > ujt p . In this figure, we find 

LO— — OC LO . (5-8) 

dio 



Therefore, the total energy, which will be obtained by 



dN 

dusLO— — , (5-9) 
dw 



is logarithmically divergent, which is consistent with the diverging radiated energy obtained in Sec. IV based on the 
point-splitting regularization. 

For u > uq, if we respect the causal structure of the background LTB spacetime, in which the naked singularity 
is ingoing null, an ingoing null ray never changes into an outgoing null ray. It is obvious that the geometrical optics 
approximation will be no longer valid after the formation of the naked singularity. Therefore, we need to determine 
mode functions by other methods than the geometrical optics approximation. On the other hand, we have seen that 
the unbounded power and energy are radiated before the occurrence of the naked singularity. Then, we could expect 
that naked singularity formation will be prevented for some physical mechanism. For example, the back reaction of 
the quantum effect might prevent naked singularity formation. In fact, the LTB spacetime is only an approximation 
of realistic gravitational collapse. Violation of spherical symmetry, rotational support, hardening of an equation of 
state, ignition of nuclear burning, quantum gravity, and so on might make the approximation by the LTB spacetime 
no longer valid in the final stage of the realistic collapse. From the above discussions, the causal structure of the LTB 
spacetime after naked singularity formation cannot be taken seriously. 



B. Momentary spectrum 

It is important in estimating the validity of the geometrical optics approximation to examine which frequency band 
dominates the power at some moment. For this purpose, we determine the momentary spectrum by wavelet analysis. 
We use Gabor's mother wavelet| in place of e ix because it has a very clear physical meaning as a Gaussian wave 
packet. We consider Gabor's mother wavelet 



''The Gabor wavelet is sometimes referred to as the Morlet wavelet. Strictly speaking, Gabor's mother wavelet does not become 
a basis function. However, it is known that this wavelet with an appropriate value of a is useful in finding the frequency of the 
signal. 
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ip(x) = 



4tt<j 



-e e 



(5.10) 



where we have set a = 8. The result below is not so sensitive to the value of a. We can obtain the Bogoliubov 
coefficients as 



1 UJ 



^y-y dvTP[u(F(v)-u)]e-™\ 
^ = -tt\ - / dvp[u>(F{v) - u)]e^'\ 



where it is chosen to be the observation time. Because the Gabor wavelet occupies area AwAu 
(u,u>), the momentary number spectrum is given by 



dN 1 
dujdu 2 



(5.11) 
(5.12) 

2 in the phase space 
(5.13) 



where r w is the transmission coefficient, which we tentatively assume — 1. 

We have calculated numerically the momentary spectrum. The result is shown in Fig. |^. This figure displays the 
contribution of each logarithmic bin of frequency to the power at the moment. The parameters are fixed as A = 0.8 
and /; = 1. We have chosen the time of observation u = IO^cj" 1 . The shape of the spectrum looks somewhat like the 
thermal one but the normalization of the spectrum is different. In this figure, it is found that the contribution to the 
power is dominated by the frequency uj ~ 2tt(uq — u) _1 . For a conformally coupled scalar field, the integration of the 
obtained ujdN / (dcodu) with u> recovers the power derived by the point-splitting regularization except for a numerical 
factor of order unity, i.e., 



dhJLO- 



dN 
(kudu 



P. 



(5.14) 



For a minimally coupled scalar field, such a correspondence cannot be seen. For a conformally coupled scalar field, 
since the power is mainly transported by particles with a frequency of about 2ir(uo — u) , the radiated number of 
particles per unit time is constant, being roughly estimated from Eq. (4.17) as 



From this equation, we can roughly estimate as 



du 

dN 
doj 



(5.15) 



(5.16) 



which recovers the total spectrum except for a numerical factor. From the above analysis, we can conclude that the 
total divergence term in the expression for the power of a minimally coupled scalar field has nothing to do with the 
observed particle number. 



VI. NULL GEODESICS AND REDSHIFTS IN LTB SPACETIME 



A. Assumption 



Here we analytically examine null geodesies and redshifts in the L TB spacetime and derive some useful res ults t o 
comprehend the problem. We assume the function F(r) of the form ( |3.8D with F3 > and F5 < 0. From Eq. (3.17), 
R and i?. r are written for r/I < 1 as 




-1/3 



(6.1) 
(6.2) 
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where r\ is defined as 

T) = t - t . 

Moreover, T) a (r) = t s {r) — to and r) a h(r) = t a h(r) — to are approximated for r/lo < 1 as 

Vs(r) 1 / r 



to 2 \IqJ 
Vahjr) _ 1 (t 
to ~ 2 



2 /r*\3/T\3/ \3 
2 \ fl \ IT 



3/ \^o/ vo 



For simplicity, we assume 



(6.3) 

(6.4) 
(6.5) 

(6.6) 



This implies that we can safely expand in powers of r/lo and take the leading-order term in the whole of the cloud. 



B. Trajectories 



First we consider the region where 



(6.7) 



is satisfied. This region is approximately recognized as the Friedmann universe. From Eq. (3.4), it is found that this 
condition can be satisfied only for r\ < 0. In this case, R^ r is approximated as 



R r 



to 



2/3 



The ordinary differential equation (4.1) can be easily integrated to 



to 



1/3 



1 fh\ (r 



3 \toJ Wo 



Ca±, 



(6.8) 



(6.9) 



where Ca± is a constant of integration. It is easily found that each null ray can be drawn by a parallel transport of 
the curve with Ca± — along the r-axis direction. 
Next we consider the region where 



lo 



(6.10) 



is satisfied. This region is not at all approximated by the Friedmann universe. In this case, R^ r is approximated as 



Then Eq. (4.1) is integrated to 



to 



R r 



±2 -2/3 



7 



2 2 / 3 3 V^o 



h 

to J \l 



4/3 



7/3 



— Cb±, 



(6.11) 



(6.12) 



where Cb± is a constant of integration. It is easily found that each null ray can be drawn by a parallel transport of 
the curve with Cb± = along the ?y-axis direction. 

See Fig. M, w hich illustrates the spacetime around [r), r) — (0,0). Condition (6.7) is satisfied in region A, while 
condition ( |b\10| ) is satisfied in region B. The boundary of regions A and B will be described by 
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to \h 



(6.13) 



where 7 is a constant of order unity. We denote this boundary curve as C. It may be kept in mind that this treatment 
is rather simple. However, we believe that the present approximation, where we divide the spacetime into two regions, 
will be enough to comprehend the essence of the problem. 

Then, we concentrate on the behavior of null geodesies around the naked singularity. See Fig. ^ for the trajectory 
of null geodesies. Fro m th e regular center 77 < and r = 0, the outgoing null ray runs region A which is described by 
the upper sign of Eq. ( |6.9| ) with Ca+ > 0. Then this null ray goes through boundary C. After that the null ray goes 
into region B and the trajectory is described by the upper sign of Eq. ( 3.12| ) with Cb+ > 0. Since all null rays with 
Cb+ > in region B emanate from the regular center, the outgoing null geodesic with C'b+ = in region B generates 
the Cauchy horizon. In other words, it is the earliest outgoing null ray from the singularity at r = 0. We denote 
this outgoing null geodesic as no+. There are infinitely many outgoing null geodesies later than 72.0+ which emanates 
from the naked singularity at r = 0. In fact, these null geodesies form one parameter family. These null rays cannot 
asymptote to no+ as r — > because these null rays are obtained by parallel transport of no+ in the ?7-axis direction 
in region B. Instead, the location of the apparent horizon r\ = rj a h(r) is an asymptote of these outgoing null rays in 
approaching the naked singularity. 

On the other hand, it is clear that the ingoing null ray with Cb- = terminates at the singularity r = 0. We 
denote this ingoing null ray as no-. For small positive Cb-, the ingoing null ray in region B crosses to boundary C 
and the null ray becomes described by the lower sign of Eq. (6i)) with Ca- > 0, and then terminates at the regular 
center. For small negative Cb-, the ingoing null ray in region B terminates at the spacelike singularity r > 0. 



C. Condition for globally naked singularity 

If no + reach 3 + , then the singularity is globally naked, otherwise the singularity is locally naked. Noting that the 
intersection of the apparent horizon with the cloud surface is on the event horizon, we find that the condition for the 
singularity to be globally naked is given by 



lo 



1/3' 



16 
27 



lo 



1/3 



(6.14) 



Since x 1 
naked as 



(16/27)a; 2 is a decreasing function of x > 0, we can find the condition for the singularity to be globally 



^<0.66(7 



(6.15) 



This implies that the singularity is globally naked for sufficiently small r s f when to and 1$ are fixed. Thus, it turns 
out that the present assumption r s f/lo <C 1 is relevant for the globally naked singularity case. We can translate the 
above condition to a condition for M and lj s as 

M < 6.4 x l(T 2 uj- 1 . (6.16) 

It implies that the mass of the dust cloud with the globally naked singularity is bounded from above. 

D. Redshifts 



Since it is necessary for later discussions, we turn our attention to redshifts of particles. Let fc M be the tangent 
vector of a radial null geodesic. Inside of the dust cloud, satisfies the following equations: 



dX 



{k l f + (Rr) 2 (k r ) 2 = 0, 
dk l 



dX 



+ R r Rtr(k r Y = 0, 



[{R, r ) 2 k r ] - i? r i? rr (fc r ) 2 = 0, 



(6.17) 
(6.18) 

(6.19) 
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where A is the affine parameter. The frequency uj observed by an observer comoving with a fluid element is calculated 
from fc M as 

uj = -k^u^ = -k l . (6.20) 

In region A, Eqs. ( |6.17 )-( |6~19 ) can be integrated as 

fc*(-ry) 2/3 ps const, (6.21) 

k r (-v) 4/3 ~ const (6.22) 

i n the lowest order. Therefore, a particle is blueshifted for the comoving observer in region A. In region B, Eqs. 
(S.19) can be integrated as 

fc* ps const, (6.23) 

k r r 4/3 ps const (6.24) 

in the lowest order. Therefore a particle is neither redshifted nor blueshifted in region B. 

In the external Schwarzschild spacetime, in a similar way, we obtain the tangent vector of a radial null geodesic as 

kT ( 1 ~ -TT ) = const > ( 6 - 25 ) 



R J 

k R = const. (6.26) 

A static distant observer observes the frequency 

uj = -fc T (T,oo). (6.27) 

Using the matching conditions ( 3.12| )-(3.14), the observed frequency uj = —k^u^ by the comoving observer at the 
surface is written using the above w as 

f=r^-, (6-28) 
1 =F a 

where a is defined by 



\ ™, (6-29) 



and R is the areal radius of the dust surface when the null ray crosses the surface. 



E. Estimate of function G or F 

Let a ± be the value of a for Uq±. Let r)a± be the value of rj when n ± crosses the dust surface. Then, from 



Eq. (3.16), it is easy to see that the Taylor expansion is valid around i] = ?/ - for v as a function of r\ as 



v = v + — (V - Vo-) + 0(( V - 77 _) 2 ) . (6.30) 

1 + «o- 

For ao+ < 1, from Eq. ( |3.15| ), it is also easy to see that the Taylor expansion is valid around 77 = t]q + for u as a 
function of r\ as 

u = u + - — (77 - 770+) + O ((77 - ?7o+) 2 ) , (6.31) 

1 - a 0+ 

or equivalently, 

Vo+ ~ V = (1 - "o+)("o - u) + 0((uq - u) 2 ). (6.32) 

See Fig. ^ again. An ingoing light ray which is close to no± enters the dust cloud in region B (a), enters region A 
(b), crosses the center (c), becomes an outgoing light ray, enters again region B (d) and leaves the dust cloud (e). In 
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order to determine G or F, we need to relate the time coordinate r\ a at which the incoming light ray enters the dust 
cloud with rj e at which the outgoing light ray leaves the dust cloud. Within the present approximation, the following 
relations hold: 



Vo- - 




-1]b 


to 




to 


=±) 


1/3 




to J 




\t0 y 


-Vc\ 


1/3 




to J 






Vo+ - 






to 







-JOb 
to 



7/6 



1/3 



1/3 



_ i 7 - 1/2 f- 



Vb 

to J \ to 



1/2 



1/2 



where r]o± are written as 



From the above relations, we find 



Then, we find for ao+ < 1 



2 -2/3 7 -7/6 



to/ V *o 



7/6 



7/3 



f?o+ - r] e ~ r?o- - ?7a- 



u = G(u) Ps v - - a ° + (uo - u), 
1 + a - 



or 



u = F(v) ps n - (wo - v). 

1 - a + 

From the above discussions, we can determine A — lim u _^ Uo G' as 

1 - a 0+ 



A 



1 + a - 



For the collapse to a black hole, from Eq. (3.15), we can express u by r\ around R — 2m at the surface as 

u ps — 4Mhi(r)h — rj) + const, 



(6.33) 
(6.34) 
(6.35) 
(6.36) 

(6.37) 
(6.38) 
(6.39) 
(6.40) 

(6.41) 
(6.42) 



where % is the time when the event horizon crosses the dust surface. On the other hand, we can easily find the 
following linear relation when an incoming light ray with the advanced time v becomes to an outgoing null ray which 
reaches the surface at the time ry as 



Vh -r/~ const x (v h - v), 



(6.43) 



where the ingoing null ray with Vh is reflected to the outgoing ray on the event horizon. Then we find Eq. ( 4.27 ) 



VII. VALIDITY OF THE GEOMETRICAL OPTICS APPROXIMATION 



The geometrical optics approximation is exact for the two-dimensional spacetime, where the line element is given 
by setting dfl 2 = in the four-dimensional line element (3.1). Therefore, the power, energy, and spectrum obtained 
in Sees. IV and |v| are exact in this sense for the two-dimensional version of the LTB spacetime. 

In the four-dimensional case, the geometrical optics approximation is only an approximation because of the existence 
of the curvature potential. The geometrical optics approximation is valid for waves with the wavelength shorter than 
the curvature radius of the spacetime geometry. In other words, the geometrical optics approximation is good if 
condition 
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a > Shrill 1 /' 



(7.1) 



is satisfied, where the hat denotes components in the local incrtial tetrad frame. The nonvanishing components of the 
Riemann tensor are as follows: 



Rfftf ~ 




Rfitf 


~ "-rift ~ 


"■tele = 






= R§iei = 




~r9§f 






RaXal 


RhXXa 


Rj,aa 





R, 



R r 



L t<f>t4> 



rd>rd 



L t<j><j>t 



L 4>tt<f> 



L 4>t<f>t 
= R 



(f>r<fir 



R,u 
R ' 
R t R tr 

RR r 



In region A, they are approximated as 



I R-ifif I 



\Rt§t$\ 



R± 
R 



\Rf§fe\ 



\Ral 



In region B, they are approximated as 



I R-ifif I 



\R 



I0I()\ 



\R. 



t 2 \ 1 



(7.2) 

(7.3) 
(7.4) 

(7.5) 
(7.6) 

(7.7) 



In Sec. VI D, we have seen that u> is kept constant approximately in region B and that il> may be considerably 
blucshiftcd in region A. However, in Sec. VIE, we can see that the null ray which goes into the dust cloud enters 
region A and get out of region A at the same time in the lowest order, i.e, r\\, w r\ c k, rjd- This implies that the 
frequency of such a particle is kept almost constant all over the dust cloud. On the other hand, the Riemann tensor 
along the light ray reaches the maximum when the light ray goes through region A. Now we can write down the 
condition for the geometric optics approximation to be valid as 



u>2t:(-t 1c )- 1 . 

This condition can be rewritten by the quantities on 3 + using Eqs. ( 6.28| ) and ( 3.31 ). The result is 

uj > u> cr (u) = 27r(uo — u)^ 1 



(7.8) 
(7.9) 



for «o+ < 1- 

As we have shown in Sec. [v|, the power for each mode is mainly transported by particles of which the frequency 
is about 2tt(uq — it)" 1 . Since the transmission coefficient T w will be a function of (to /uj cr ), it is natural to estimate 
r w (w = uj cr ) as of the order unity. It implies that the calculations based on the geometrical optics approximation will 
be valid for a rough order estimate. 



VIII. DISCUSSIONS 



Here we compare the naked singularity explosion with the famous radiative phenomena, the Hawking radiation, 
and the black hole evaporation. The Hawking radiation is derived in the fixed background spacetime as has been 
done in this paper. The result is thermal radiation, the temperature of which is given by 



and the power is given by 



^= 768^' (8 - 2) 



The Hawking radiation is not explosive but constant. 

In contrast to the naked singularity explosion, it is not until the back reaction of the quantum effect is taken into 
account that explosive radiation is emitted. For the black hole evaporation, the power is given by 
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p~(* e „-t)- 2/3 , 



(8.3) 



where t ev is the time of evaporation and the above equation is valid only for t ev — t > 1, where 1 means, of course, 
the Planck time. The divergence is apparently much milder than the naked singularity explosion. The expression of 
the power does not have any parameter dependence. The radiated energy during the final evaporation phase is only 
finite in contrast to the naked singularity. The mass of the evaporating black hole is given by 



M ~ (U-t) 1 / 3 . 

From the above equations, the power at the time t is dominated by particles with frequency 

LO ~ (t ev - t)- 1 ^. 



The total spectrum is estimated from Eqs. (3.3) and 

dN 

LO 

duo 



5ft as 



(8.4) 



•5) 



(8.6) 



where this estimate is valid for lo < 1. Equation (3.6) implies that the spectrum is apparently softer than that of the 
naked singularity explosion. For the Hawking radiation, since the Schwarzschild black hole has only one parameter, 
i.e., its mass, the back reaction effect is rather easy to be taken into account although it is phenomenological. For the 
naked singularity, the system is dynamical and not parametrized by only one parameter. It suggests that it is not an 
easy task to take the back reaction of the quantum effect into account even phenomenologically. The comparison of 
radiation processes in the collapse of a dust ball is shown in Table ffl. 

TABLE I. Comarison of radiation processes: Hawking radiation (HR), black hole evaporation (BHE), naked singularity 
explosion (NSE) Q, and naked singularity explosion in self-similar spacetime (NSESS) |lq-|2C 



Process 


HR 


BHE 


NSE 


NSESS 


Radiation 


Constant 


Explosive 


Explosive 


Explosive 


P (total div. subtracted) 


(7687TM 2 )- 1 


^ (tev £) ^ 


~ Vs(t C H - t)^ 1 


~ (icH — t)~ 2 


E 




M 


~ V a log(*CH - t) 


~ (tcH — i)" 1 


dN 
dv 


V 






? 


e 167r 2 Mi/ ^ 






< v > 


~ M -1 


dM_ 


~ (£ch — t) 


? 


Back reaction 




dt 
Good 






Geometrical approx. 


Good 


Marginal 


? 


Redshift 


oo 


oo 


Finite 


oo 
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Again we return to the naked singularity explosion. The singularity frequency io s is obtained from Eq. ( |4.9| ) as 

^ „ io-3 ( AS) "7 *2_ V Hz . (8.7) 
Vims/ V 10km / 

The characteristic frequency depends not only on the free fall time of the system but also greatly on the scale 
of inhomogeneity. The scale of inhomogeneity becomes larger, the singularity frequency becomes higher. In some 
situations, gravitational collapse will set in due to the growth of small perturbations of hydrostatic equilibrium in a 
Newtonian regime. In such a case, to and lo will be related as lo/to = \/l5c s , where c s is the sound speed at the 
center. Then the singularity frequency is given by uj s ~ 10 3 c^tg 1 . It implies that lj s does not much exceed the inverse 
of the initial free fall time in such situations. 



From Eq. (|5.16|) , it is found that 

JuHp duJ U tp 



Thus, the total number of emitted particles with frequency above wt p is a few at most. It suggests that the radiation 
may be extremely directional because a few particles transport almost all the energy radiated during the collapse. 

We are not sure how to take into account the back reaction of the quantum effect. However, we can assume that 
the back reaction of the quantum effect becomes important if a considerable fraction of the total energy of the system 
are radiated away by the quantum particle creation. Fro m this assumption, we can estimate the emitted energy as 



E ~ 1O 54 (M/M0)erg, although M must satisfy condition (6T6). For E — fM with / < 1, the averaged energy u a 
of the radiated particles above iot v is given by 

toav ~ E/N ~ E{Mu a )- 2 '* ~ fM x ' 3 u)- 2 ' 3 . (8.9) 
The central energy density e c at the final stage is given by 

e c ~ u% ex P (—) ~ M-^loT exp (—) , (8.10) 



because At is given by 



Us J \u s J 



Af ~<exp(-^) ~A/ 2 /V 1/2 exp(-^') . (8 . n) 



Actually, this estimate for the radiated energy should be considered as an upper bound. It is noted that we have 
to transform these intrinsic quantities to observed quantities by considering possible reactions of created energetic 
particles. 

The naked singularity explosion may be a new candidate for a source of ultrahigh energy cosmic rays because it 
has a potentially extremely hard spectrum. It may also provide a candidate for a central engine of 7-ray bursts. The 
naked singularity explosion could contribute 7-ray background radiation or disturb big bang nucleosynthesis. 



IX. CONCLUSIONS 



We have considered the naked singularity which occurs in the collapse of a dust ball with the initial density 
profile which is physically reasonable. We have calculated radiation during the naked singularity formation under 
the geometrical optics approximation. Since the regularization in these calculations has been done only in a flat 
spacetime, the results contain no ambiguity which could come from local curvature. We have obtained the power and 
energy of the radiation, which turn out to diverge to positive infinity for both minimally and conformally coupled 
massless scalar fields. We have also obtained the spectrum of radiation, both total and momentary ones. Then we 
have found that the power is dominantly contributed from particles of frequency to ~ 27r(icH — t) , where (ten — t) is 
the remaining time until a distant observer could see the singularity. The spectrum is given by iodN/du) « (gj/u) s ) , 
which is much harder than that of the black hole evaporation. The radiation depends on the characteristic frequency 
lu s of the singularity, which is determined by the initial free fall time and the initial scale of inhomogeneity at the 
center. The radiation may be expected to be strongly beamed. The geometrical optics approximation has turned out 
to be consistent for a rough order estimate. It will provide a very interesting candidate for a central engine of 7-ray 
bursts or a source of ultrahigh energy cosmic rays because the predicted intrinsic spectrum may be extended to an 
extremely high frequency. A detailed study on such astrophysical implications will be a future work. 

After this manuscript was received, the paper j34j appeared, in which the naked singularity formation is discussed 
as a possible origin of 7-ray bursts. 
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APPENDIX A: CONSERVED QUANTITY IN LTB COLLAPSE 



In Sec. [II, we set the initial data at t — and specified the solution. In fact, the choice of an initial slice is arbitrary. 
Therefore, another slice among the synchronous comoving slicings can be chosen as an initial slice. Since this variety 
of choices of an initial slice, of course, does not change the original collapse model, physical quantities, such as, power 
and energy, must not depend on the choic e of an in itial slice. 

At an arbitrary time t < to, using Eqs. (3.8) and (6.1), we can expand the mass function F(r) in terms of the areal 
radius R as 



I?- 5 



-13/3 



(Al) 



We can define the coefficients of R 3 and R 5 as -^3(77) and (77) . Then we can define the free fall time £0(77) and the 
scale of inhomogeneity Iq(j]) at the time 77 as 



f (7,)^^( J F 3 (7,))- 1 / 2 =io 
-1/2 



to 



h(v) 



-h(v) 



hi-? 

to 



7/6 



(A2) 
(A3) 



where to and are understood as to(0) and in(0), respectively. Then we can easily find that the following quantity 
is independent of the choice of an initial slice. 



to(v) 7 to 
= -77 = const. 

^6 It 



(A4) 



APPENDIX B: ANALYTIC EXPRESSIONS OF BOGOLIUBOV COEFFICIENTS 



As seen in Sec. VI, for ao+ < 1, we can write the function F(v) in the following form 

F(v) « uo - a(vo -v) + b(v - v) s , 

where 5, a and b satisfy 

a > 0, 
b > 0. 

We denote the advanced time as v = vo—d after which the higher-order terms can be neglected, where d - 
is satisfied. Because the particle creation will occur due to the third term in the right-ha nd sid e of 
dominant contribution to the spectrum comes from the integr al on Vp — d < v < v in Eqs. ( 2.15f ) and 
the Bogoliubov coefficients are obtained from Eqs. ( 2.15 ) and ( 2.16 ) as 



(Bl) 




27T V it) 

1 fu 



g—iiiJUQ—iu> f VQ I ^iuj(as — bs )-\-iu)' ' s 



(B5) 
(B6) 
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They have the following expression 
1 



'*VVf-l#-^8*— — — / dxe 

2^2^ L )* k \ ° ( w ')*+(*-i)(M L axe 



27T\/ LOU} 1 

1 



k=0 j=0 

oc k 

-iuuQ — iuj' vq 



c k 



Cl> C 



Z^Z^ ^ fc! - ( ^k+is-i^k-j) j 

k=0j=0 ' 



27TV OJOj' 

where we have introduced the frequency c which characterizes the singularity as 

c=bW- l \ 

They also have another expansion: 
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They also have another expansion: 
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For lu' — ► 0, as seen in Eqs. (B7) and (B8), they behave as 



27TV WW' 

i 



For a/ — > oo, as seen in Eqs. ([BIO] ) and ( pll|) , they behave as 

1 (w') 1/2 
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where we have used the formula 



{e + ia) z 



(B7) 
(B8) 

(B9) 



(BIO) 



(Bll) 

(B12) 
(B13) 

(B14) 
(B15) 

(B16) 
(B17) 

(B18) 



for e > and z > 0. 
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FIG. 1. Null rays interior of the dust cloud in the LTB spacetime. 
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FIG. 2. Plots of (a) G(u), and (b) G'(u), (c) G"(w) for the LTB spacetime. 
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FIG. 3. Power for minimally and conformally coupled scalar fields in the LTB spacetime 
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FIG. 4. Null rays interior of the dust cloud in the Oppenheimer-Snyder spacetime. 
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FIG. 5. Plots of (a) G(u), (b) G'(it), and (c) G"(u) for the Oppenheimer-Snyder spacetime. 
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FIG. 7. Total spectrum of a naked singularity explosion for A = 0.8 and /; = 1. 
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FIG. 8. Momentary spectrum of a naked singularity explosion at tto — u — lO -6 ^ 1 for A = 0.8 and /, = 1. 
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FIG. 9. Schematic diagram of the LTB spacetime interior of the dust cloud around the naked singularity 
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